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A data sample corresponding to an integrated luminosity of 2.1 fb _1 collected by the DO detector 
at the Fermilab Tevatron Collider was analyzed to search for squarks and gluinos produced in pp 
collisions at a center-of-mass energy of 1.96 TeV. No evidence for the production of such particles 
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was observed in topologies involving jets and missing transverse energy, and 95% C.L. lower limits of 
379 GeV and 308 GeV were set on the squark and gluino masses, respectively, within the framework 
of minimal supergravity with tan/3 = 3, Ao = 0, and /i < 0. The corresponding previous limits are 
improved by 54 GeV and 67 GeV. 

PACS numbers: 14.80.Ly, 12.60.Jv, 13.85.Rm 



In supersymmetric models eacn °f the standard 
model (SM) particles has a partner differing by a half- 
unit of spin. If i?-parity Q is conserved, supersymmetric 
particles are produced in pairs, and their decay leads to 
SM particles and to the lightest supersymmetric particle 
(LSP), which is stable. Cosmological arguments suggest 
that the LSP should be neutral and colorless [3(. The 
lightest neutralino Xii which is a mixture of the super- 
partners of the neutral gauge and Higgs bosons, fulfills 
these conditions. In the following, it is assumed that im- 
parity is conserved and that x? is the LSP. Since this 
particle is weakly interacting, it escapes detection and 
provides the classic missing transverse energy ($t) signa- 
ture at colliders. In pp collisions, squarks (q) and gluinos 
(g), the superpartners of quarks and gluons, would be 
abundantly produced, if sufficiently light, by the strong 
interaction, leading to final states containing jets and 
1£t- The most constraining direct limits on squark and 
gluino masses were obtained by the DO collaboration [J] , 
based on an analysis of 310 pb -1 of data from pp col- 
lisions at a center-of-mass energy of 1.96 TeV, collected 
during Run II of the Fermilab Tevatron Collider. In the 
model of minimal supergravity (mSUGRA) [5] , the limits 
obtained were rriq > 325 GeV and m~ g > 241 GeV at the 
95% C.L., for the set of model parameters detailed below. 
In this Letter, an update of the DO search for squarks and 
gluinos in topologies with jets and large J^t is reported, 
using a seven times larger data set of 2.1 fb _1 . 

A detailed description of the DO detector can be found 
in Ref. Q . The central tracking system consists of a sil- 
icon microstrip tracker and a central fiber tracker, both 
located within a 2 T superconducting solenoidal magnet. 
A liquid-argon and uranium calorimeter covers pseudora- 
pidities up to \rj\ k 4.2, where r\ = — In [tan (6/2)], and 9 
is the polar angle with respect to the proton beam direc- 
tion. The calorimeter consists of three sections, housed 
in separate cryostats: the central one covers \r]\ < 1.1, 
and the two end sections extend the coverage to larger 
\r)\. The calorimeter is segmented in depth, with four 
electromagnetic layers followed by up to five hadronic 
layers. It is also segmented in semi-pro jective towers of 
size 0.1 x 0.1 in the (77, <fi) plane, where 4> is the azimuthal 
angle in radians. Calorimeter cells are defined by the 
intersections of towers and layers. Additional sampling 
is provided by scintillating tiles between cryostats. An 
outer muon system, covering \rj\ < 2, consists of a layer 
of tracking detectors and scintillation trigger counters in 
front of 1.8 T iron toroids, followed by two similar lay- 
ers after the toroids. Jets were reconstructed with the 



iterative midpoint cone algorithm [7j with cone radius 
1Z = y/ (A0) 2 + (Ay) 2 = 0.5 in azimuthal angle <\> and 
rapidity y = \ \n((E + p z )/{E - p z )). The jet energy 
scale (JES) was derived from the transverse momentum 
balance in photon-plus-jet events. The J^t was calcu- 
lated from all calorimeter cells, and corrected for the jet 
energy scale and for the transverse momenta of recon- 
structed muons. 

The DO trigger system consists of three levels, LI, 
L2 and L3. The 1.2 ftr 1 of data recorded after 2006 
(during the so-called Run lib) were collected using a 
significantly upgraded system, in particular a new LI 
calorimeter trigger Q involving a sliding-window algo- 
rithm which improved the jet triggering efficiency. In ad- 
dition, JfiT was used at LI to select events, which was not 
done during the previous data taking period (Run Ha). 
Both improvements helped to keep a high trigger effi- 
ciency despite the increased instantaneous luminosity in 
Run lib (up to 2.8 x 10 32 cm _2 s~ 1 , to be compared to 
1.6 x 10 32 cm~ 2 s _1 in Run Ha). The events used in this 
analysis were recorded using two categories of triggers 0] . 
The dijet triggers selected events with two jets and IpT, 
while the multijet triggers were optimized for events with 
at least three jets and 

The SM processes leading to events with jets and 
real Ifix in the final state ("SM backgrounds") are the 
production of W or Z bosons in association with jets 
(W/Z+jets), of pairs of vector bosons (WW, WZ, ZZ) 
or top quarks (tt), and of single top quarks. The neu- 
trinos from the decays Z — ► vv and W — + Iv, with the 
W boson produced directly or coming from a top quark 
decay, generate the Jfir signature. In this analysis, most 
of the W boson leptonic decays leading to an electron 
or a muon were identified, and the corresponding events 
rejected. However, a charged lepton from W boson de- 
cay can be a tau decaying hadronically. It can also be 
an electron or a muon that escapes detection or fails the 
identification criteria. Such M^+jets events therefore ex- 
hibit the jets plus 1£t signature. Finally, multijet pro- 
duction also leads to a final state with jets and when 
one or more jets are mismeasured ("QCD background"). 

Events from SM and supersymmetric processes were 
simulated using Monte Carlo (MC) generators and passed 
through a full GEANT3-based [l0( simulation of the de- 
tector geometry and response. They were subsequently 
processed with the same reconstruction chain as the data. 
The parton density functions (PDFs) used in the MC 
generators are the CTEQ6L1 [IH PDFs. A data event 
from a randomly selected beam crossing was overlaid on 
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each event to simulate the additional minimum bias in- 
teractions. The QCD background was not simulated, but 
estimated directly from data. The ALPGEN generator fl2| 
was used to simulate W/Z+jets and tt production. It 
was interfaced with PYTHIA [l3} for the simulation of 
initial and final state radiation (ISR/FSR) and of jet 
hadronization. Pairs of vector bosons and electroweak 
top quark production were simulated with pythia and 
COMPHEP [14j, respectively. The next-to-leading order 
(NLO) cross sections were computed with mcfm 5.1 [HI]. 

Squark and gluino production and decay were simu- 
lated with pythia. The masses of the supersymmetric 
particles were calculated with SUSPECT 2.3 [l6| from the 
set of five mSUGRA parameters: mo, the universal scalar 
mass at the grand unification scale Aqut; fn\i% the uni- 
versal gaugino mass at Aqut; Aq, the universal trilinear 
coupling at Aqut; tan/3, the ratio of the vacuum ex- 
pectation values of the two Higgs fields; and \i the sign 
of the Higgs-mixing mass parameter. The decay widths 
and branching ratios of all supersymmetric particles were 
then calculated with SDECAY 1.1a 17]. In order to allow 
for an easier comparison with previous results, the fol- 
lowing parameters were fixed to the same values as in 
Ref. [J]: A) = 0, tan/3 = 3, and fi < 0. The produc- 
tion of scalar top quarks (stops) was ignored, and in the 
following, the "squark mass" is the average mass of all 
squarks other than stops. All squark and gluino decay 
modes were taken into account, including cascade decays 
through charginos or neutralinos which could addition- 
ally produce electrons, muons, or taus. The NLO cross 
sections of the squark and gluino pair production were 
calculated with prospino 2 jlgj . 

The analysis strategy is the same as in Ref. [3] with 
three analyses optimized for three benchmark regions 
of the mSUGRA parameter space. A "dijet" analysis 
was optimized at low mo for events containing a pair of 
acoplanar jets, as expected from pp — > qq — > qXiQXi an( i 
pp — > qq — > QXiQXv A "gluino" analysis was optimized 
at high mo for events with at least four jets, as expected 
from pp — > gg — > qqXiQQXv Finally, a "3-jets" analysis 
was optimized for events with at least three jets, as ex- 
pected from pp — > qg — > qXi<l<jXi- The benchmark for 
this analysis is the case where m§ = m~ g . 

Each analysis required at least two jets and substan- 
tial IpT (> 40 GeV). The acoplanarity, i.e. the azimuthal 
angle between the two leading jets, was required to be 
smaller than 165°, where the two leading jets are those 
with the largest transverse momenta. These leading 
jets were required: to be in the central region of the 
calorimeter, 1 77det | < 0.8, where r^ot is the jet pseudo- 
rapidity calculated under the assumption that the jet 
originates from the detector center; to have fractions of 
energy in the electromagnetic layers of the calorimeter 
less than 0.95; and to have transverse momenta greater 
than 35 GeV. The best primary vertex (PV0) was se- 
lected among all reconstructed primary vertices (PV) 



TABLE I: Selection criteria for the three analyses (all energies 
and momenta in GeV); see the text for further details. 



Preselection Cut 




All Analyses 








> 40 




Vertex z pos.[ 




< 60 cm 




Acoplanarity 




< 165° 




Selection Cut 


"dijet" 


"3-jets" 


"gluino" 


Trigger 


dijet 


multijet 


multijet 


jeti pr a 


> 35 


> 35 


> 35 


jet 2 p T a 


> 35 


> 35 


> 35 


j et 3 PT b 


— 


> 35 


> 35 


ipf. . t)it> 

J Cb 4 PI 






> 20 


Electron veto 


yes 


yes 


yes 


Muon veto 


yes 


yes 


yes 


A0(#r,jeti) 


> 90° 


> 90° 


> 90° 


A</>(# T ,jet 2 ) 


> 50° 


> 50° 


> 50° 


A(frmin{$T, anyjet) 


> 40° 






Ht 


> 325 


> 375 


> 400 




> 225 


> 175 


> 100 



"First and second jets are also required to be central (|)7d e t I < 0.8), 
with an electromagnetic fraction below 0.95, and to have 
CPF0 > 0.75. 

''Third and fourth jets are required to have l^jetl < 2.5, with an 
electromagnetic fraction below 0.95. 



as the one with the smallest probability to be due to 
a minimum bias interaction |19(. Its longitudinal posi- 
tion with respect to the detector center was restricted to 
\z\ < 60cm to ensure efficient reconstruction. For each 
jet, CPF0 is defined as the fraction of track pr sum as- 
sociated with the jet which comes from PV0, CPF0= 
£>)p ck (PV0)/X>T ack ( an y PV). The two leading jets 
were required to have CPF0 larger than 0.75. Since jet 
transverse momenta and Ipx were calculated with respect 
to PV0, this criterion reduced the background with large 
Ijjf due to an incorrect PV0 selection. 

Different selection criteria were applied in the three 
analyses, as summarized in Table HI Events passing a 
trigger for acoplanar dijet events were used in the "dijet" 
analysis, while a multijet and $t trigger was required in 
the "3-jets" and "gluino" analyses. In the "3-jets" and 
"gluino" analyses, a third and fourth jet were required, 
respectively. In comparison with the previously pub- 
lished analysis the |?7dct| upper limit for the third and 
fourth jet was increased from 0.8 to 2.5, which consider- 
ably improved the signal efficiency at high mo. These two 
jets also had to have their fraction of energy in the elec- 
tromagnetic layers of the calorimeter less than 0.95, but 
no requirement was made on their CPF0. The third and 
fourth jet transverse momenta were required to exceed 
35 GeV and 20 GeV, respectively. In all three analyses, 
a veto on isolated electrons or muons with pt > 10 GeV 
was applied to reject background events containing a lep- 
tonic W boson decay. The azimuthal angles between the 
$t and the first jet, A0($r,jeti), and the second jet, 
A0(^?t, jet2), were used to remove events where the en- 
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ergy of one jet was mismeasured, generating ]£t aligned 
with that jet. In the "dijet" analysis, the minimum az- 
imuthal angle A0 m j n (^r, anyjet) between the ]£t and 
any jet with p? > 15 GeV was required to be greater 
than 40° to further suppress QCD background. This cri- 
terion was not used in the "3-jets" and "gluino" analyses 
because of the higher jet multiplicity in the signal events. 

The two final cuts on ]£t and on Ht — SjetsPT, where 
the sum is over all jets with px > 15 GeV and \rjdet I < 2.5, 
were optimized by minimizing the expected upper limit 
on the cross section in the absence of signal. Here and in 
the calculation of the final limits, the modified frequentist 
CL S method [2(| was used. The QCD background contri- 
bution was estimated by fitting the $t distribution below 
60 GeV with an exponential function, after subtraction of 
the SM background processes, and subsequently extrap- 
olating this function above the chosen ]£t cut value. The 
optimal cuts are given in Table HI The fx distributions 
after applying all analysis criteria except the one on IpT 
are shown in Fig. [T] for the three analyses. 

Table HI] reports the number of data events and the 
expected signal and background. The main background 
contributions are from [Z — > vv) + jets, (W — > Iv) +jets, 
and tt-^bbqq'lv. The QCD background was evaluated 
from a fit to the $t distribution as described above. The 
largest QCD contribution of 1.4±0.8 event was estimated 
in the "gluino" analysis, but was conservatively ignored 
when setting the limits. It was found to be negligible in 
the "dijet" and "3-jets" analyses. The signal efficiencies 
are given in Table HT1 for the three benchmark scenarios, 
with the corresponding values of mo, mi/2, the squark 
and gluino masses, and the NLO cross section. 

The uncertainty coming from the JES corrections is 
typically (10-15)% for the SM backgrounds and (6-11)% 
for the signal efficiencies. The uncertainties due to the 
jet energy resolution, to the jet track confirmation, and 
to jet reconstruction and identification efficiencies range 
between 2% and 4%. All these uncertainties on jet prop- 
erties account for differences between data and MC sim- 
ulation, both for signal efficiencies and background con- 
tributions. The trigger was found to be fully efficient 
for the event samples surviving all analysis cuts with an 
uncertainty of 2%. The uncertainty on the luminosity 
measurement is 6.1% 21] . All of these uncertainties are 



fully correlated between signal and SM backgrounds. A 
15% systematic uncertainty was set on the W/Z+jets 
and ti NLO cross sections. The uncertainty on the signal 
acceptance due to the PDF choice was determined to be 
6%, using the forty-eigenvector basis of the CTEQ6.1M 
PDF set [III]. Finally, the effects of ISR/FSR on the 
signal efficiencies were studied by varying the pythia 
parameters controlling the QCD scales and the maximal 
allowed virtualities used in the simulation of the space- 
like and time-like parton showers. The uncertainty on 
the signal efficiencies was determined to be 6%. 

The nominal NLO signal cross sections, a nom , were 



computed with the CTEQ6.1M PDF and for the renor- 
malization and factorization scale p r j = Q, where Q was 
taken to be equal to nig for gg production, rriq for qq and 
qq production, and (rriq + mg)/2 for qg production. The 
uncertainty due to the choice of PDF was determined 
using the full set of CTEQ6.1M eigenvectors, with the 
individual uncertainties added in quadrature. The ef- 
fect on the nominal signal cross sections, which varies 
between 15% and 60%, is dominated by the large uncer- 
tainty on the gluon distribution at high x. The effect 
of the renormalization and factorization scale was stud- 
ied by calculating the signal cross sections for /x r ,f = Q, 
Mr,f = Q/2 and p, r ,f = 2 x Q. The factor of two on 
this scale reduces or increases the nominal signal cross 
sections by (15-20)%. The PDF and n r j effects were 
added in quadrature to compute minimum, a nim , and 
maximum, er ma x, signal cross sections. 

The numbers of events observed in the data are in 
agreement with the SM background expectation in the 
three analyses. Therefore, an excluded domain in the 
gluino-squark mass plane was determined as follows. The 
three analyses were run over signal MC samples gener- 
ated in the gluino-squark mass plane to compute signal 
efficiencies. Then, to take advantage of the different 
features of the three analyses, seven independent com- 
binations of the three selections were defined as shown 
in Table III, which were combined in the limit compu- 
tations. The number of events passing at least one of 
the three analyses is 31 while the SM expectation is 
32.6 ± 1.7(stat.)±^(syst.) events. Figure M shows the 
95% C.L. observed and expected upper limits on squark- 
gluino production cross sections for the three benchmark 
scenarios. Figure [3] shows the excluded domain in the 
gluino-squark mass plane. 

The absolute lower limits on the squark and gluino 
masses obtained in the most conservative hypothesis, 
cr min , are 379 GeV and 308 GeV, respectively. The corre- 
sponding expected limits are 377 GeV and 312 GeV. Ta- 
ble [IV] summarizes these absolute limits for the three sig- 
nal cross section hypotheses. Limits were also derived for 
the particular case rriq = m g . For cr m i n , squark and gluino 
masses below 390 GeV are excluded, while the expected 
limit is 390 GeV. The observed limit becomes 408 GeV 
for cr nom , and 427 GeV for cr max . 

The results of this analysis also constrain the 
mSUGRA parameters at the grand unification scale. Fig- 
ure|4]shows the excluded regions in the (mo, m^) plane 
for tan/? — 3, A = 0, /i < 0. Although a detailed scan 
of the mSUGRA parameter space is beyond the scope of 
this analysis, it was verified that similar results hold for 
a large class of parameter sets. In particular, the fact 
that there is no explicit veto in this analysis on hadron- 
ically decaying tau leptons mitigates the expected effi- 
ciency reduction at larger values of tan /3. It can be seen 
in Fig. Q] that the limits from LEP2 chargino and slepton 
searches [23j | are improved for tuq values between 70 and 
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FIG. 1: Distributions of 1$t after applying all analysis criteria except the one on IpT for the "dijet" (a), "3-jets" (b), and "gluino" 
(c) analyses; for data (points with error bars), for SM background (full histogram), and for signal MC (dotted histogram on top 
of SM). The signal drawn corresponds to (mo,mi/2) = (25, 175) GeV, rriq — m§ = 410 GeV, and (mo, rai/2) = (500, 110) GeV 
from left to right. The fitted QCD background contribution is also shown. 



TABLE II: For each analysis, information on the signal for which it was optimized (mo, myj, m g , m,-, and nominal NLO cross 
section), signal efficiency, the number of events observed, the number of events expected from SM backgrounds, the number 
of events expected from signal, and the 95% C.L. signal cross section upper limit. The first uncertainty is statistical and the 
second is systematic. 

Analysis (m ,m 1/2 ) (m g ,mq) <r nom e sig . AW AWkgrd. A^ig. 095 

(GeV) (GeV) (pb) (%) (pb) 

"dijet" (25,175) (439,396) 0.072 6.8 ±0.4^;;; 11 11.1 ± 1-2±|;1 10.4±0.6±^ 0.075 

"3-jets" (197,154) (400,400) 0.083 6.8 ±0.4^;* 9 lO.7iO.9l2 J 12.0 ± 0.7±|;| 0.065 

"gluino" (500,110) (320,551) 0.195 4.1 ±0.3^ j 20 17.7 ± l-l±f ;| 17.0 ± 1.2±|;| 0.165 



TABLE III: Definition of the analysis combinations, and number of events observed in the data and expected from the SM 
backgrounds. 



Selection 


"dijet" 


"3-jets" 


"gluino" 


AW 


-^backgrd. 


Combination 1 


yes 


no 


no 


8 


9.4 ± 1.2 (stat.) ±\% (syst.) 


Combination 2 


no 


yes 


no 


2 


4.5 ± 0.6 (stat.) (syst.) 


Combination 3 


no 


no 


yes 


14 


12.5 ± 0.9 (stat.) t\% (syst.) 


Combination 4 


yes 


yes 


no 


1 


1.1 ± 0.3 (stat.) ±g"| (syst.) 


Combination 5 


yes 


no 


yes 




kinematically not allowed 


Combination 6 


no 


yes 


yes 


4 


4.5 ± 0.6 (stat.) t\% (syst.) 


Combination 7 


yes 


yes 


yes 


2 


0.6 ± 0.2 (stat.) ±° l (syst.) 


At least one selection 








31 


32.6 ± 1.7 (stat.) (syst.) 




FIG. 2: For tan/3 = 3, Ao = 0, fj, < 0, observed (closed circles) and expected (opened triangles) 95% C.L. upper limits on 
squark-gluino production cross sections combining the analyses for mo = 25 GeV (a), rriq = m g (b), and mo = 500 GeV (c). 
The nominal production cross sections are also shown, with shaded bands corresponding to the PDF and renormalization-and- 
factorization scale uncertainties. 
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TABLE IV: Absolute lower limits at the 95% C.L. on the 
squark and gluino masses (in GeV) as a function of the choice 
of signal cross section hypothesis as defined in the text. Num- 
bers in parentheses correspond to the expected limits. These 
limits are valid for the mSUGRA parameters tan/3 = 3, 
A a = 0, /x < 0. 



Hypothesis 


Gluino mass 


Squark mass 


^min 


308 (312) 


379 (377) 


^nom 


327 (332) 


392 (391) 


^max 


349 (354) 


406 (404) 



LEP2 f 



D0, L=2.1 fb" 
tanp=3, A q =0, |^<0 




no EWSB 



100 200 300 400 500 600 

m (GeV) 




Gluino Mass (GeV) 

FIG. 3: In the gluino and squark mass plane, excluded regions 
at the 95% C.L. by direct searches in the mSUGRA frame- 
work with tan/3 = 3, Ao = 0, < 0. The region excluded by 
this analysis and previous DO Run II results [4} in the most 
conservative hypothesis (er m i n ) is shown in dark shading. The 
thick (dotted) line is the limit of the observed (expected) ex- 
cluded region for the cr n om hypothesis. The band delimited by 
the two dashed lines shows the effect of the PDF choice and of 
the variation of /x r ,f by a factor of two. Regions excluded by 
previous experiments are indicated in light shading (2^. The 
two thin lines indicate the indirect limits inferred from the 
LEP2 chargino and slepton searches The region where 
no mSUGRA solution can be found is shown hatched. 



300 GeV and for miw values between 125 and 165 GeV. 
However, the LEP2 Higgs search limits remain more con- 
straining in a purely mSUGRA scenario [231 ] . 

In summary, a search for squarks and gluinos pro- 
duced in pp collisions at 1.96 TeV has been performed in 
a 2.1 fb _1 data sample. The results of three selections of 
events with jets and large missing transverse energy are 
in agreement with the SM background predictions. In 
the framework of minimal supergravity with tan/3 = 3, 



FIG. 4: In the (mo, TO1/2) plane, the region excluded by this 
analysis at the 95% C.L. in the mSUGRA framework for 
tan/3 = 3, Ao = 0, /i < is shown in dark shading. The 
thick line is the limit of the excluded region for the cr noln hy- 
pothesis. The corresponding expected limit is the dashed line. 
The band delimited by the two dotted lines shows the effect of 
the PDF choice and of the variation of /x r ,f by a factor of two. 
Regions excluded by the LEP2 chargino and slepton searches 
are indicated in light shading [23| ]. The region where there is 
no electroweak symmetry breaking is shown in black. 

A Q = 0, and fj, < 0, 95% C.L. lower limits of 392 GeV and 
327 GeV were set on the squark and gluino masses, re- 
spectively, for the central choice of PDF and for a renor- 
malization and factorization scale equal to the mass of 
the squark or gluino produced. Taking into account the 
PDF uncertainties and allowing for a factor of two in 
the choice of scale, these limits are reduced to 379 GeV 
and 308 GeV. They exceed the corresponding previous 
limits [4[ by 54 GeV and 67 GeV and are the most con- 
straining direct limits on the squark and gluino masses 
to date. 
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